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We describe two experimental tests of the Equivalence Principle that are based on frequency 
measurements between precision oscillators and/or highly accurate atomic frequency standards. 
Based on comparisons between the hyperfine frequencies of 87 Rb and 133 Cs in atomic fountains, 
the first experiment constrains the stability of fundamental constants. The second experiment is 
based on a comparison between a cryogenic sapphire oscillator and a hydrogen maser. It tests Local 
Lorentz Invariance. In both cases, we report recent results which improve significantly over previous 
experiments. 

PACS numbers: 



I. INTRODUCTION 



Einstein's Equivalence Principle (EEP) is at the heart of special and general relativity and a cornerstone of modern 
physics. The central importance of this postulate in modern theory has motivated tremendous work to experimentally 
test EEP 0. Additionally, nearly all unification theories (in particular string theories) violate EEP at some level 
0- El which further motivates experimental searches for such violations. A third motivation comes from a recent 
analysis of absorption systems in the spectra of distant quasars p( which seems to indicate a variation of the fine- 
structure constant a over cosmological timescale, in violation of EEP. 

EEP equivalence principle is made of three constituent elements. The Weak Equivalence Principle (WEP) postu- 
lates that trajectories of neutral freely falling bodies are independent of their structure and composition. Local Lorentz 
Invariance (LLI) postulates that in any local freely falling reference frame, the result of a non gravitational measure- 
ment is independent of the velocity of the frame. Finally, Local Position Invariance (LPI) postulates that in any local 
freely falling reference frame, the result of a non gravitational measurement is independent of where and when it is 
performed. The experiments described here use precision oscillators and atomic clocks to test LLI and LPI. 

In a first section, experiments testing LPI are described. In these experiments, frequencies of atomic transitions 
are compared to each other in a local atomic clock comparison. The measurements are repeated over a few years. 
LPI implies that these measurements give consistently the same answer, a prediction which is directly tested. These 
experiments can be further interpreted as testing the stability of fundamental constants, if one recognizes that any 
atomic transition frequency can (at least in principle) be expressed as a function of properties of elementary particles 
and parameters of fundamental interactions. Such interpretation requires additional input from theoretical calculations 
of atomic frequencies. In this first section, after reviewing some of these calculations, we describe a comparison between 
87 Rb and 133 Cs hyperfine frequencies in atomic fountains. We give new and significant constraints to the stability of 
fundamental constants based on the results of these experiments. 

In a second part, a test of LLI is described. The frequency of a macroscopic cryogenic sapphire resonator is 
measured against a hydrogen maser as a function of time. Wc look for sidereal and and semi-sidereal modulations 
of the measured frequency which would indicate a violation of LLI depending on the speed and orientation of the 
laboratory frame with respect to a preferred frame. First the Robertson, Mansouri and Sexl theoretical framework is 
described as a basis to interpret the experiments. Then new results improving on previous experiments are reported. 
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II. TEST OF LOCAL POSITION INVARIANCE. STABILITY OF FUNDAMENTAL CONSTANTS 



A. Theory 

Tests described here are based on highly accurate comparisons of atomic energies. In principle at least, it is 
possible to express any atomic energy as a function of the elementary particle properties and the coupling constants 
of fundamental interactions using Quantum Electro-Dynamics (QED) and Quantum Chromo-Dynamics (QCD). As a 
consequence, it is possible to deduced a constrain to the variation of fundamental constants from a measurement of 
the stability of the ratio between various atomic frequencies. 

Different types of atomic transitions are linked to different fundamental constants. The Iryperfine frequency in a 
given electronic state of alkali-like atoms (involved for instance in 133 Cs, 87 Rb Q, 199 Hg + [21 13, 171 Yb + microwave 
clocks) can be approximated by: 

^RooCxMxg^(^) a*F^(a), (1) 



'hfs - ±L oo^ ~ ~ h f s ~ y y— j « ^ h fs 

where Roc is the Rydberg constant, c the speed of light, grW the nuclear g-factor, m e /m p the electron to proton mass 
ratio and a the fine structure constant. In this equation, the dimension is given by RooC, the atomic unit of frequency. 
-^Ms i s a numerical factor which depends on each particular atom. Fj^ s (a) is a relativistic correction factor to the 
motion of the valence electron in the vicinity of the nucleus. This factor strongly depends on the atomic number 
Z and has a major contribution for heavy nuclei. The superscript (i) indicates that the quan tity depends on each 
particular atom. Similarly, the frequency of an electronic transition (involved in H 9] , 40 Ca P3, 199 Hg+ [HI, 171 Yb+ 
|12l [T3 optical clocks) can be approximated by 

V<L* Boo* * JUL (2) 

Again, the dimension is given by RooC A^ cc is a numerical factor. F^ c (a) is a function of a which accounts for 
relativistic effects, spin-orbit couplings and many-body effects [65|. From equations ^ and |3 it is possible to calculate 
the ratio between the frequencies in atomic species (i) and (ii), depending on the type of transition involved: 

Hi) 

^clcc „ ^clccW ( n\ 
'clcc r clec\ a ) 

^O.g^^^sM (4) 

„('*) Mi) F {ll) (a\ 

(5) 

The product of constants RoqC cancels out in the ratio of two atomic frequencies and only dimensionless factors are 
left. Also, numerical factors that are irrelevant to the present discussion have been omitted. Already, the different 
sensitivity of the various type of comparison can be seen in these equations. Comparisons between electronic transitions 
(equation |3J are only sensitive to a. Comparisons between hyperfine transitions (equation [SJ are sensitive both to 
a and the strong interaction through the nuclear g-factors. Comparisons between an electronic transition and a 
hyperfine transition (equation^) are sensitive to a, to the strong interaction and to the electron mass. 

The sensitivity of a given atomic transition to the variation of fundamental constants can be derived from equation 
[2 and 12 
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In recent work [Ijj, it has been suggested that this analysis can be pushed one step further by linking the g-factors 
giW and the proton mass m p to fundamental parameters, namely the mass scale of QCD Aqcd and the quark mass 
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Tfiq — (m u + m<j)/2. Within this framework, any atomic frequency comparison can be interpreted as testing the 
stability of 3 dimensionless fundamental constants: a, m q j Aqcd and m e /AQcrj- For any transition we can write: 



S l n ( ^L) * jfW x ^ + K W x ^Mqcd) + ^) x ^K/Aqcd) (g) 



-RooC/ "a 9 (m g /A QC D) ' (m e /A QC D) 

where the superscript (z) again indicates that the coefficient depends on each particular transition. Hyperfine tran- 
sitions are sensitive to all three fundamental constants (Ka,Kq ^ 0; ~ 1). For electronic transitions, we 

have Ke,Kq ~ and therefore they are essentially sensitive to variation of a. Four well chosen atomic transitions 
constraining the stability of 3 independent frequency ratios are enough to constrain independently the stability of 
the three fundamental constants a, to 9 /Aqcd and m e /AQCD- From these equations, it can be seen that at least two 
different hyperfine transitions must be involved to independently constrain m„/ Aor; n and m e /Aq)CD which empha- 
sizes the need to maintain and improve highly accurate microwave atomic clocks 66] . With more than 4 well chosen 
atomic clocks, redundancy is achieved which means that a non vanishing variation of fundamental constants can be 
identified by a clear signature. 

Calculations of the coefficients have now been done for a large number of atomic species 0, 0, ITHl IT^. UtI llH IT^ | . 
For hyperfine transitions in 87 Rb, 133 Cs and 199 Hg + the most recent calculation gives [l4j: 

§hx / *W 87 Rb)\ „ 2M 5a _ a064 ^K/A Q cD) + 5(m e /A QCB ) 

V RoaC J a (m g /AQCDj (me/AQCDj 

5ln ( ^r° s)> ) * 2.83^ + 0.11 5 P4P°4 + (10) 
V -RooC / a (m 9 /A QCD ) (TO e /A QCD ) 

51n / ^ fs ( 199 Hg + ) \ ^ i 3 S_a _ ^(m g /A QCD ) J(m e /A QCD ) 



V -RcoC / a (m g /AQCD) (™ e /AQCD) 

In each case^the most recent and precise value for the K$ coefficient given here is in good agreement with earlier 
calculations y|- For electronic transitions in H, 40 Ca and 199 Hg + , we have: 

ikf^l- 0. (12) 
§ln (^fC & )\ 0Q3 Sa (13) 



RoaC J a 

8ln[^ r HS+) V-3-2^ d4) 



RocC ) a 

A new generation of laser cooled optical clocks is now under development in several groups ( 87 Sr |13,|^,|22|, 171 Yb, 
27 A1 + [23j,...). This work will significantly improve the stringency and redundancy of this test of LPI. 



B. Experiments with Rb and Cs fountain clocks 

In these experiments, three atomic fountains are compared to each other, using a hydrogen maser (H-maser) as a 
flywheel oscillator (Fig^J. Two fountains, a transportable fountain FOM, and FOl |2J| are using cesium atoms. The 
third fountain is a dual fountain (F02) [2j|, operating alternately with rubidium (F02Rt>) and cesium (F02c s ). These 
fountains have been continuously upgraded in order to improve their accuracy from 2 x 10~ 15 in 1998 to 8 x 10~ 16 
for cesium and from 1.3 x 10~ 14 |26j to 6 x 10~ 16 for rubidium. 

Fountain clocks operate as follows. First, atoms are collected and laser cooled in an optical molasses or in a 
magneto-optical trap in 0.3 to 0.6 s. Atoms are launched upwards, and selected in the clock level (tuf — 0) by 
a combination of microwave and laser pulses. Then, atoms interact twice with a microwave field tuned near the 
hyperfine frequency, in a Ramsey interrogation scheme. The microwave field is synthesized from a low phase noise 
100 MHz signal from a quartz oscillator, which is phase locked to the reference signal of the H-maser (FigUJ. After 
the microwave interactions, the population of each hyperfine state is measured using light induced fluorescence. This 
provides a measurement of the transition probability as a function of microwave detuning. Successive measurements 
are used to steer the average microwave field to the frequency of the atomic resonance using a digital servo system. 
The output of the servo provides a direct measurement of the frequency difference between the H-maser and the 
fountain clock. 
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FIG. 1: BNM-SYRTE clock ensemble. A single 100 MHz signal from a H-maser is used for frequency comparisons and is 
distributed to each of the microwave synthesizers of the 133 Cs (FOl, FOM, F02 Cs ) and 87 Rb fountain clocks. In 2001, the Rb 
fountain has been upgraded and is now a dual fountain using alternately rubidium (F02Rb) or cesium atoms (F02c s ). 



The three fountains have different geometries and operating conditions: the number of detected atoms ranges 
from 3 x 10 5 to 2 x 10 6 at a temperature of ~ 1 /iK, the fountain cycle duration from 1.1 to 1.6 s. The Ramsey 
resonance width is between 0.9 and 1.2 Hz. In measurements reported here the fractional frequency instability is 
(1 — 2) x 10 _13 t -1 / 2 , where t is the averaging time in seconds. Fountain comparisons have a typical resolution of 
~ 10~ 15 for a 12 hour integration, and each of the four data campaigns lasts from 1 to 2 months during which an 
accuracy evaluation of each fountain is performed. 




52560 52570 52580 52590 52600 
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FIG. 2: The 2002 frequency comparison data, a) H-maser fractional frequency offset versus FOM (black squares), and al- 
ternately versus F02Rb (open circles) and F02c s (open triangles between dotted lines), b) Fractional frequency differences. 
Between dotted lines, Cs-Cs comparisons, outside Rb-Cs comparisons. Error bars are purely statistical. They correspond to the 
Allan standard deviation of the comparisons and do not include contributions from fluctuations of systematic shifts of Table 

The 2002 measurements are presented in Fig|21 which displays the maser fractional frequency offset, measured by 
the Cs fountains FOM and F02c s - Also shown is the H-maser frequency offset measured by the Rb fountain F02u,b 
where the Rb hyperfme frequency is conventionally chosen to be fRb(1999) = 6 834 682 610.904 333 Hz, our 1999 value. 
The data are corrected for the systematic frequency shifts listed in Table [I] The H-maser frequency exhibits fractional 
fluctuations on the order of 10~ 14 over a few days, ten times larger than the typical statistical uncertainty resulting 
from the instability of the fountain clocks. In order to reject the H-maser frequency fluctuations, the fountain data 
are recorded simultaneously (within a few minutes). The fractional frequency differences plotted in Fig|5]b illustrate 
the efficiency of this rejection. F02 is operated alternately with Rb and Cs, allowing both Rb-Cs comparisons and 
Cs-Cs comparisons (central part of Fig|2J) to be performed. 

Systematic effects shifting the frequency of the fountain standards are listed in Table [I] The quantization magnetic 
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TABLE I: Accuracy budget of the fountains involved in the 2002 measurements (F02 et FOM). 



Fountain 


F02 Cs 


F02 Rb 


FOM 


Effect 


Value & Uncertainty (10~ ib ) 


2 nd order Zeeman 


1773.0 ± 5.2 


3207.0 ±4.7 


385.0 ±2.9 


Blackbody Radiation 


-173.0 ±2.3 


-127.0 ± 2.1 


-186.0 ±2.5 


Cold collisions 
+ cavity pulling 


-95.0 ±4.6 


0.0 ± 1.0 


-24.0 ±4.8 


others 


0.0 ±3.0 


0.0 ±3.0 


0.0 ±3.7 


Total uncertainty 


8 


6 
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field in the interrogation region is determined with a 0.1 nT uncertainty by measuring the frequency of a linear field- 
dependent "Zeeman" transition. The temperature in the interrogation region is monitored with 5 platinum resistors 
and the uncertainty on the black-body radiation frequency shift corresponds to temperature fluctuations of about 1 
K 1271 • Clock frequencies are corrected for the cold collision and cavity pulling frequency shifts using several methods 
|28l l29j|. For Rb, unlike [2j|, an optical molasses with a small number of atoms (~ 5.4 x 10 6 ) is used. We thus 
estimate that these two shifts are smaller than 5 x 1CP 17 . All other effects do not contribute significantly and their 
uncertainties are added quadratically. We searched for the influence of synchronous perturbations by changing the 
timing sequence and the atom launch height. To search for possible microwave leakage, we changed the power (x9) in 
the interrogation microwave cavity. No shift was found at a resolution of 10 -15 . The shift due to residual coherences 
and populations in neighboring Zeeman states is estimated to be less than 10 -16 . As shown in jsij, the shift due to 
the microwave photon recoil is very similar for Cs and Rb and smaller than +1.4 x 10~ 16 . Relativistic corrections 
(gravitational redshift and second order Doppler effect) contribute to less than 10 -16 in the clock comparisons. 
For the Cs-Cs 2002 comparison, we find: 

^(2002) -^"(2002) = x iQ i6 



fCs 

where the first parenthesis reflects the la statistical uncertainty, and the second the systematic uncertainty, obtained 
by adding quadratically the inaccuracies of the two Cs clocks (see Table HJ. The two Cs fountains are in good 
agreement despite their significantly different operating conditions (see Table [1}, showing that systematic effects are 
well understood at the 10 -15 level. 

In 2002, the 87 Rb frequency measured with respect to the average 133 Cs frequency is found to be: 

2/R b (2002) = 6 834 682 610.904 324(4) (7) Hz (16) 

where the error bars now include F02ptb, F02c s and FOM uncertainties. This is the most accurate frequency 
measurement to date. 

In FigElare plotted all our Rb-Cs frequency comparisons. Except for the less precise 1998 data pflj . two Cs fountains 
were used together to perform the Rb measurements. The uncertainties for the 1999 and 2000 measurements were 
2.7 x 10~ 15 , because of lower clock accuracy and lack of rigorous simultaneity in the earlier frequency comparisons 
|25l |. A weighted linear fit to the data in FigOU determines how our measurements constrain a possible time variation 
of i^Rb/^Cs- We find: 

d .l n (^BL] =(0.2 ±7.0) xlO- 1 ^ 1 (17) 



dt \ fc, 

which represents a 5-fold improvement over our previous results pfij and a 100-fold improvement over the Hg + -H 
hyperfine energy comparison Q ■ Using equation [5] and the sensitivity to a in equation |T5] and ^3 we find that this 
results implies the following constraint: 

(/ • ( 9Cs^ „pA9\ _ (n o -i- 7 a\ v m-i6, rT .-i 



■In =(0.2 + 7.0) x lO-^yr- 1 . (18) 

dt \gvLb 



Using the link between g- factors, m q and Aqcd (ref. .14j, eq. 19*1 and ITU)l . we get: 

d 



■In 
dt 



(a 49 [m g /A QCD ] 017 ) = (0.2 ± 7.0) x Hr 18 yr _1 - (19) 
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FIG. 3: Measured 87 Rb frequencies referenced to the 133 Cs fountains over 57 months. The 1999 measurement value (i/Rb(1999) = 
6 834 682 610.904 333 Hz) is conventionally used as reference. A weighted linear fit to the data gives j- t In = (0.2 ± 7.0) x 

10~ 16 yr -1 . Dotted lines correspond to the la slope uncertainty. 

A comparison between the single mercury 199 Hg + ion optical clock and the 133 Cs hyperfine splitting has been 
recently reported by the NIST group [llj which (according to eq. ^| and I14|) constrains the stability of 
a 6 [m e /A Q cD][m g /A Q cD] 01 at the level of 7 x 1CT 15 yr" 1 . 



III. TESTS OF LOCAL LORENTZ INVARIANCE 



Local Lorentz Invariance is one of the constituent elements of the Einstein Equivalence Principle (see Sect. 1) and 
therefore one of the cornerstones of modern physics. It is the fundamental hypothesis of special relativity and is related 
to the " constancy of the speed of light" . The central importance of this postulate in modern physics has motivated 
tremendous work to experimentally test LLI 0, l3l| . Additionally, nearly all unification theories (in particular string 
theory) violate the EEP at some level M which further motivates experimental searches for such violations of the 
universality of free fall Q and of Lorentz invariance [33l l34| . 

We report here on experimental tests of LLI using a cryogenic sapphire oscillator and a hydrogen maser. The 
relative frequency of the two clocks is monitored looking for a Lorentz violating signal which would modulate that 
frequency at, typically, sidereal and semi-sidereal periods due to the movement of the lab with the rotation of the 
Earth. We set limits on parameters that describe such Lorentz violating effects, improving our previous results [3jJ 
by a factor 2 and the best other results [3^, ^| by up to a factor 70. 

Many modern ex peri ments that test LLI rely essentially on the stability of atomic clocks and macroscopic resonators 
[35L l36t I37I l38l 13^. l4Ct l4l| , therefore improvements in oscillator technology have gone hand in hand with improved 
tests of LLI. Our experiment is no exception, the large improvements being a direct result of the excellent stability of 
our cryogenic sapphire oscillator. Additionally its operation at a microwave frequency allows a direct comparison to 
a hydrogen maser which provides a highly stable and reliable reference frequency. 

Numerous test theories that allow the modeling and interpretation of experiments that test LLI have been developed. 
Kinematical frameworks |42l |43| postulate a simple parametrisation of the Lorentz transformations with experiments 
setting limits on the deviation of those parameters from their special relativistic values. A more fundamental approach 
is offered by theories that parametrise the coupling between gravitational and non-gravitational fields (THe/j, [ll l44ll45| 
or xg 46] formalisms) which allow the c omp arison of experiments that test different aspects of the EEP. Finally, 
formalisms based on string theory I31L I3^l33| have the advantage of being well motivated by theories of physics 
that are at present the only candidates for a unification of gravity and the other fundamental forces of nature. 



A. Theory 

Owing to their simplicity the kinematical frameworks o f l42l |43| have been widely used to model and interpret 
many previous experiments testing LLI j3f|, |37|, |H, |3^, 0, l48| and we will follow that route. An analysis based on 
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the more fundamental "Standard Model Extension" (SME) is under way and will be published shortly. 

By construction kinematical frameworks do not allow for any dynamical effects on the measurement apparatus. 
This implies that in all inertial frames two clocks of different nature (e.g. based on different atomic species) run at 
the same relative rate and two length standards made of different materials keep their relative lengths. Coordinates 
are defined by the clocks and length standards, and only the transformations between those coordinate systems are 
modified. In general this leads to observable effects on light propagation in moving frames but, by definition, to no 
observable effects on clocks and length standards. In particular, no attempt is made at explaining the underlying 
physics (e.g. modified Maxwell and/or Dirac equations) that could lead to Lorentz violating light propagation but 
leave e.g. atomic energy levels unchanged. On the other hand dynamical frameworks (e.g. the THe/i formalism or the 
SME) in general use a modified general Lagrangian that leads to modified Maxwell and Dirac equations and hence to 
Lorentz violating light propagation and atomic properties, which is why they are considered more fundamental and 
more complete than the kinematical frameworks. Furthermore, as shown in [49|. the SME is kept sufficiently general 
to, in fact, encompass the kinematical frameworks and some other dynamical frameworks (in particular the THe/i 
formalism) as special cases, although there are no simple and direct relationships between the respective parameters. 

Concerning our experiment the SME allows the calculation of Lorentz violating effects on the fields inside the 
sapphire resonator, on the properties of the sapphire crystal itself and on the hydrogen maser transition. As shown 
in [59j the effect on the sapphire crystal amounts to only a few percent of the direct effect on the fields, and [60j show 
that the hydrogen mp — — > m' F = clock transition is not affected to first order. Hence the total effect is dominated 
by the Lorentz violating properties of the electromagnetic fields inside the resonator which can be calculated following 
the principles laid down in 49] . An analysis of our experiment in that framework is currently being carried out, and 
the results will be published elsewhere in the near future. In this paper we concentrate on the analysis using the 
kinematical framework of Mansouri and Sexl 431 . 



The Robertson, Mansouri &: Sexl Framework 

Kinematical frameworks for the description of Lorentz violation have been pioneered by Robertson [42^ ] and further 
refined by Mansouri and Sexl ^3| and others. Fundamentally the different versions of these frameworks are equivalent, 
and relations between their parameters are readily obtained. As mentioned above these frameworks postulate gener- 
alized transformations between a preferred frame candidate £(T, X) and a moving frame S(i,x) where it is assumed 
that in both frames coordinates are realized by identical standards. We start from the transformations of (in 
differential form) for the case where the velocity of S as measured in £ is along the positive X-axis, and assuming 
Einstein synchronization in S (we will be concerned with signal travel times around closed loops so the choice of 
synchronization convention can play no role): 

,™ 1 ( , vdx\ dx v ( , vdx\ dy ,_ dz . . 

dT=-[dt + ^);dX = — + -[dt + —^)-,dY=-r;dZ=— 20 
a \ cr J b a \ c z ) d d 

with c the velocity of light in vacuum in E. Using the usual expansion of the three parameters (a ~ 1 + cims 1 * 2 /c 2 + 
0(4); b sa 1 + f3 MS v 2 /c 2 + 0(4); d « I + 5 MS v 2 /c 2 + 0(4)), setting c 2 dT 2 = dX 2 + dY 2 + dZ 2 in S, and transforming 
according to l|20[l we find the coordinate travel time of a light signal in S: 

dl { v 2 (\ \ v 2 \ 

dt=- M -(ftjg-OMs-l)— - {--f3MS + 5 M sUm 2 e—\ +0(4) (21) 



where dl = y 7 dx 2 + dy 2 + dz 2 and 6 is the angle between the direction of light propagation and the velocity v of S 
in S. In special relativity cims = —1/2; /3ms = 1/2; 5ms = an d <|2*0|l reduces to the usual Lorentz transformations. 
Generally, the best candidate for £ is taken to be the frame of the cosmic microwave background (CMB) [5(1151] with 
the velocity of the solar system in that frame taken as u© ~ 377 km/s, deck rj —6.4°, RA k, 11. 2h. 

Michelson-Morley type experiments H^, [5^ determine the coefficient Pmm — (1/2 — /?ms + 5ms) of the direction 
dependent term. For many years the most stringent limit on that parameter was \Pmm\ < 5 x 10~ 9 determined 
over 23 years ago in an outstanding experiment |35|. Our experiment confirms that result with roughly equivalent 
uncertainty (2.2 x 10~ 9 ). Recently an improvement to |Pmm| < 1-5 x 10~ 9 has been reported |4l|. Kennedy-Thorndikc 
experiments 1361 l37L l38j measure the coefficient Pkt = (Pms — &ms — 1) of the velocity dependent term. The most 
stringent limit |38[]on |Prtt| has been recently improved from [3^ by a factor 3 to |Pxr| < 2.1 x 10~ 5 . We improve this 
result by a factor of 70 to |Pxt| < 3.0 x 10~ 7 . Finally clock comparison and Doppler experiments 0,E3j[l3 measure 
QfMSi currently limiting it to |q;ms + 1/2 1 < 8 x 10~ 7 . The three types of experiments taken together then completely 
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TABLE II: Present limits on Lorentz violating parameters in the framework of |43| 



Reference a M s + 1/2 1/2- ,3ms + 5ms /3ms - «ms - 1 

[4£, 4£,£2] < 8 x 1(T 7 

[25J - < 5 x 10~ 9 

[41] - (2.2 ± 1.5) x 10" 9 

[22] - - (1.9 ±2.1) x 1CT 5 

our previous results [39] - (1.5 ± 4.2) x 1CT 9 (-3.1 ± 6.9) x 1CT 7 

this work - (1.2 ± 2.2) x 10 -9 (1.6 ± 3.0) x 1(T 7 



characterize any deviation from Lorentz invariance in this particular test theory, with present limits summarized in 
Table |nl 

Our cryogenic oscillator consists of a sapphir e crystal of cylindrical shape operating in a whispering gallery mode 
(see Fig. 0]for a schematic drawing and |541.l55| for a detailed description). Its coordinate frequency can be expressed 
by v = m/t c where t c is the coordinate travel time of a light signal around the circumference of the cylinder (of radius 
r) and m is a constant. From (|21[1 the relative frequency difference between the sapphire oscillator and the hydrogen 
maser (which, by definition, realizes coordinate time in S |56|) is 

= Pkt — — + Pmm — —7T sin 9{t, <p)dtp + 0(3) (22) 

Vq C C ZTT Jq 

where vq = m/(2Trr/c), v(t) is the (time dependent) speed of the lab in E, and tp is the azimuthal angle of the 
light signal in the plane of the cylinder. The periodic time dependence of v and 9 due to the rotation and orbital 
motion of the Earth with respect to the CMB frame allow us to set limits on the two parameters in 1)22(1 by adjusting 
the periodic terms of appropriate frequency and phase (see [Hsj for calculations of similar effects for several types of 
oscillator modes). Given the limited durations of our data sets (< 15 days) the dominant periodic terms arise from 
the Earth's rotation, so retaining only those we have v(t) = u + Jl x R with u the velocity of the solar system with 
respect to the CMB, £1 the angular velocity of the Earth, and R the geocentric position of the lab. We then find after 




FIG. 4: Typical relative frequency stability of the CSO - H-maser difference after removal of a linear frequency drift. The inset 
is a schematic drawing of the cylindrical sapphire oscillator with the Poynting vector P in the whispering gallery (WG) mode, 
the velocity v(t) of the cylinder with respect to the CMB, and the relevant angles for a photon in the WG mode. 
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some calculation. 



Ai//V = P K T(HsinX) . . 

+Pmm (AcosA + 5cos(2A) + CsinA + Z?sinAcosA + £sinAcos(2A)) ^ ' 

where A = fit + <fi, and A-E and <fi are constants depending on the latitude and longitude of the lab (« 48.7°N and 
2.33°E for Paris). Numerically H « -2.6 x 10~ 9 , A « -8.8 x 1CT 8 , B « 1.8 X 1CT 7 , C-E of order 1CT 9 . We note 
that in ill'-ll) the dominant time variations of the two combinations of parameters are in quadrature and at twice the 
frequency which indicates that they should decorelate well in the data analysis allowing a simultaneous determination 
of the two (as confirmed by the correlation coefficients given in Sect. 3.2). Adjusting this simplified model to our 
data we obtain results that differ by less than 10% from the results presented in Sect. 3.2 that were obtained using 
the complete model f(|22|l including the orbital motion of the Earth). 



B. Experimental Results 

The cryogenic sapphire oscillator (CSO) is compared to a commercial (Datum Inc.) active hydrogen maser whose 
frequency is also regularly compared to caesium and rubidium atomic fountain clocks in the laboratory |25j. The 
CSO resonant frequency at 11.932 GHz is compared to the 100 MHz output of the hydrogen maser. The maser signal 
is multiplied up to 12 GHz of which the CSO signal is subtracted. The remaining w 67 MHz signal is mixed to 
a synthesizer signal at the same frequency and the low frequency beat at « 64 Hz is counted, giving access to the 
frequency difference between the maser and the CSO. The instability of the comparison chain has been measured and 
does not exceed a few parts in 10 16 . The typical stability of the measured CSO - maser frequency after removal of a 
linear frequency drift is shown in Fig. Since September 2002 we are taking continuous temperature measurements 
on top of the CSO dewar and behind the electronics rack. Starting January 2003 we have implemented an active 
temperature control of the CSO room and changed some of the electronics. As a result the diurnal and semi-diurnal 
temperature variations during measurement runs (s=s 2 weeks) were greatly reduced to less than 0.025° C in amplitude 
(best case), and longer and more reliable data sets became available. 

Our previously published results are based on data sets taken between Nov. 2001 and Sep. 2002 which did not 
all include regular temperature monitoring and control. Here we use only data that was permanently temperature 
controlled, 13 data sets in total spanning Sept. 2002 to Aug. 2003, of differing lengths (5 to 16 days, 140 days in 
total). The sampling time for all data sets was 100 s except two data sets with r = 12 s. To make the data more 
manageable we first average all points to To = 2500 s. For the data analysis we simultaneously adjust an offset and a 
rate (natural frequency drift, typically » 1.7 x 10 -18 s _1 ) per data set and the two parameters of the model l|22|) . In 
the model (|22f) we take into account the rotation of the Earth and the Earth's orbital motion, the latter contributing 
little as any constant or linear terms over the durations of the individual data sets are absorbed by the adjusted offsets 
and rates. 

When carrying out an ordinary least squares (OLS) adjustment we note that the residuals have a significantly 
non-white behavior as one would expect from the slope of the Allan deviation of Fig. The power spectral density 
(PSD) of the residuals when fitted with a power law model of the form S y (f) = kf ^ y ields typically [i » —1.5 to 
—2. In the presence of non-white noise OLS is not the optimal regression method [6l],|62] as it can lead to significant 
underestimation of the parameter uncertainties |6l| . 

An alternative method is weighted least squares (WLS) [63] which allows one to account for non-white noise 
processes in the original data by pre- multiplying both sides of the design equation (our equation lt22|) plus the offsets 
and rates) by a weighting matrix containing off diagonal elements. To determine these off diagonal terms we first 
carry out OLS and adjust the S y (f) = kf^ power law model to the PSD of the post- fit residuals determining a 
value of /i for each data set. We then use these u values to construct a weighting matrix following the method of 
fractional differencing described, for example, in [61j. Figure shows the resulting values of the two parameters (Pkt 
and Pmm) for each individual data set. A global WLS fit of the two parameters and the 13 offsets and drifts yields 
Pmm = (1-2 ± 1.9) x 10 -9 and Pkt = (1-6 ± 2.3) x 10~ 7 (la uncertainties), with the correlation coefficient between 
the two parameters less than 0.01 and all other correlation coefficients < 0.06. The distribution of the 13 individual 
values around the ones obtained from the global fit is well compatible with a normal distribution (\ 2 = 10.7 and 14.6 
for Pmm and Pkt respectively). 

Systematic effects at diurnal or semi-diurnal frequencies with the appropriate phase could mask a putative sidereal 
signal. The statistical uncertainties of Pmm and Pkt obtained from the WLS fit above correspond to sidereal and 
semi-sidereal terms (from Q23|l) of » 7 x 10~ 16 and w 4 x 10 -16 respectively so any systematic effects exceeding these 
limits need to be taken into account in the final uncertainty. We expect the main contributions to such effects to 
arise from temperature, pressure and magnetic field variations that would affect the hydrogen maser, the CSO and 
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the associated electronics, and from tilt variations of the CSO which are known to affect its frequency. Measurements 
of the tilt variations of the CSO show amplitudes of 4.6 firad and 1.6 /irad at diurnal and semi-diurnal frequencies. 

To estimate the tilt sensitivity we have intentionally tilted the oscillator by « 5 mrad off its average position which 
led to relative frequency variations of w 3 x 10~ 13 from which we deduce a tilt sensitivity of « 6 x 10 _17 /xrad _1 . 
This value corresponds to a worst case scenario as we expect a quadratic rather than linear frequency variation for 
small tilts around the vertical. Even with this pessimistic estimate diurnal and semi-diurnal frequency variations due 
to tilt do not exceed 3 x 10~ 16 and 1 x 10~ 16 respectively and are therefore negligible with respect to the statistical 
uncertainties. 

In December 2002 we implemented an active temperature stabilization inside an isolated volume (« 15m 3 ) that 
included the CSO and all the associated electronics. The temperature was measured continously in two fixed locations 
(behind the electronics rack and on top of the dewar). For the best data sets the measured temperature variations 
did not exceed 0.02/0.01 °C in amplitude for the diurnal and semi-diurnal components. In the worst cases (the 
two 2002 data sets and some data sets taken during a partial air conditioning failure) the measured temperature 
variations could reach 0.26/0.08 °C. When intentionally heating and cooling the CSO lab by « 3°C we see frequency 
variations of « 5 x 10~ 15 per °C. This is also confirmed when we induce a large sinusoidal temperature variation 
(w 1.5°C amplitude). Using this we can calculate a value for temperature induced frequency variations at diurnal 
and semi-diurnal frequencies for each data set, obtaining values that range from w 1.3 x 10~ 15 to « 5 x 10~ 17 . 

The hydrogen maser is kept in a dedicated, environmentally controlled clock room. Measurements of magnetic field, 
temperature and atmospheric pressure in that room and the maser sensitivities as specified by the manufacturer allow 
us to exclude any systematic effects on the maser frequency that would exceed the statistical uncertainties above and 
the systematic uncertainties from temperature variations in the CSO lab. 

Our final uncertainties (the error bars in Fig. |SJ are the quadratic sums of the statistical uncertainties from 
the WLS adjustment for each data set and the systematic uncertainties calculated for each data set from l|23|) and 
the measured temperature variations. For the global adjustment we average the systematic uncertainties from the 
individual data sets obtaining ±1.2 x 10 -9 on Pmm and ±1.9 x 10 -7 on Pkt- Adding these quadratically to the 
WLS statistical uncertainties of the global adjustment we obtain as our final result Pmm = (1.2 ± 2.2) x 10~ 9 and 
Pkt = (1-6 ± 3.0) x 10~ 7 (la uncertainties). 



IV. CONCLUSION AND OUTLOOK 



We have reported on two different tests of the Einstein Equivalence Principle (EEP) using the comparison of atomic 
clocks with different atomic species on one hand, and the comparison of an atomic clock and a cryogenic sapphire 
cavity oscillator on the other. The two experiments are interpreted as testing Local Position Invariance (LPI) and 
Local Lorentz Invariance (LLI) respectively which are both constituent elements of the EEP. 
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FIG. 5: Values of the two parameters (Pkt and Pmm) from a fit to each individual data set (blue diamonds) and a global fit 
to all the data (red squares). For comparison our previously published results 39] are also shown (green triangles). The error 
bars indicate the combined uncertainties from statistics and systematic effects. 
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The test of LPI is based on the comparison of the hyperfine transitions in 87 Rb and 133 Cs using atomic fountains 
that presently reach uncertainties of (6 — 8) x 10~ 16 . Such measurements were repeated over the last 5 years to 
search for a time variation that would indicate a violation of LPI. Our present results limit a linear variation to 

3J In O^^j = (0.2 ± 7.0) x 10 _16 yr _1 which represents a 5-fold improvement over our previous results and a 

100- fold improvement over the Hg + -H hyperfine energy comparison 0|. When interpreting the results as a limit on 
the time variation of fundamental constants (c.f. Sect. 2.1) we obtain 



0.49^ + 0.17 TO * A 



a m qA 



<7xl0~ lb yr (24) 



where m q \ stands for m g /AQCD- By itself this experiment limits the time variation of a combination of two of the 
three fundamental constants of Sect. 2.1. The 199 Hg + to 133 Cs comparisons by the NIST group provide 



6.0- + 0.1— 9 — + — 



a m q A m e A 



<7xlO" 1J yr" i (25) 



where m e \ stands for m e /AQCD- Combining these two results we have two constraints on the variation of the three 
fundamental constants, with the missing third constraint requiring the comparison over time with a fourth atomic 
transition (c.f. Sect. 2.1). 

The test of LLI is based on the comparison of a hydrogen maser clock to a cryogenic sapphire microwave cavity. 
This experiment simultaneously constrains two combinations of the three parameters of the Mansouri and Sexl test 
theory (previously measured individually by Michelson-Morley and Kennedy-Thorndike experiments). We obtain 
5ms — Pms + 1/2 = 1.2(1.9) (1.2) x 1CF 9 which is of the same order as the best previous results H^EI, and /3ms — 
«ms — 1 = 1-6(2. 3)(1. 9) x 10~ 7 which improves the best previous limit by a factor of 70 (the first bracket indicates 
the lcr uncertainty from statistics the second from systematic effects). We improve our own previous results j3^ 
by about a factor 2 due to more and longer data sets and to improved temperature control of the experiment. We 
note that our value on (5ms — Pms + 1/2 is compatible with the slightly significant recent result of [41| who obtained 
<5ms - /?ms + 1/2 = (2.2 ± 1.5) x 10" 9 . 

As a result of our experiment the Lorentz transformations are confirmed in this particular test theory with an 
overall uncertainty of < 8 x 10~ 7 limited now by the determination of qms from Doppler and clock comparison 
experiments [421 H9 • This is likely to be improved in the coming years by experiments such as ACES (Atomic Clock 
Ensemble in Space 63]) that will compare ground clocks to clocks on the international space station aiming at a 10 
fold improvement on the determination of «ms- 

In the future, the two tests of LPI and LLI presented here are expected to further improve due to improvements in 
the accuracies of the atomic clocks involved and due to new experimental strategies, ultimately leading to space-borne 
versions of the experiments. 

Ongoing efforts are expected to improve the accuracy of both 87 Rb and 133 Cs to the 10 -16 level. The corresponding 
limit to variations of fundamental constants will then be decreased to ~ 10 -16 yr^ 1 or less. Recent advances in the field 
of optical frequency metrology will probably lead optical frequency standards to surpass microwave clocks. Comparing 
such standards to each other will provide very stringent limits to the variation of the fine structure constant a. To 
keep the constraints to the variation of TOq/Aqcd and m e /AQCD at the same level, further efforts and new methods 
will have to be invented to improve microwave clocks. These tests will also greatly benefit from a new generation of 
time/frequency transfer at the 10 -16 level which is currently under development for the ESA space mission ACES 
which will fly ultra-stable clocks on board the international space station in 2006 J63| and a similar project conducted 
by NASA. These missions will allow highly precise comparisons between clocks developed in distant laboratories and 
based on different atomic species and/or different technologies. 

Concerning the test of LLI, new proposals have been made to use two orthogonal resonators or two orthogonal 
modes in the same sapphire resonator placed on a rotating platform |58| . Such a set-up is likely to improve the tests 
of LLI by several orders of magnitude as the relevant time variations will now be at the rotation frequency (« 0.01 — 0.1 
Hz) which is the range in which such resonators are the most stable (w 100 fold better stability). Additionally many 
systematic effects are likely to cancel between the two orthogonal oscillators or modes and the remaining ones are 
likely to be less coupled to the rotation frequency than to the sidereal frequencies used in our experiment. Ultimately, 
it has been proposed [64| to conduct these tests on board an Earth orbiting satellite, again with a potential gain of 
several orders of magnitudes over current limits. 
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